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^ ' Abstract 

' In supersymmetric models with Dirac neutrino masses, a left-right mixed sneutrino 

can be a viable dark matter candidate. We examine the MSSM+Pr parameter space 

, where this is the case with particular emphasis on light sneutrinos with masses below 

CN ' 10 GeV. We discuss implications for direct and indirect dark matter searches, including 

the relevant uncertainties, as well as consequences for collider phenomenology. 

& 1 Introduction 

A simple extension of the Standard Model (SM) by right-handed neutrinos provides the 
framework for describing neutrino masses and the observed neutrino oscillations [H[2] (for 
O \ reviews see, e.g., [2111] )• Current observations, however, do not allow to establish the Majo- 
rana or Dirac nature of neutrinos. While the smallness of the neutrino mass can be naturally 
O ■ explained by introducing Majorana mass terms and making use of the see-saw mechanism, 
tyQ . Dirac masses for neutrinos with very small Yukawa couplings are a viable and interesting 
alternative. In supersymmetric models, one may naturally obtain very light Dirac neu- 
trino masses from F-term SUSY breaking [5]. In addition to providing an explanation for 
neutrino masses, this class of supersymmetric models offers an interesting alternative dark 
matter (DM) candidate, the sneutrino. Indeed in these models one can generate a weak- 
scale trilinear A„ term that is not proportional to the small neutrino Yukawa couplings. Thus 
large mixing between left-handed (LH) and right-handed (RH) sneutrinos can be induced 
even though the Yukawa couplings are extremely small. This is in sharp contrast with the 
usual MSSM where the trilinear A terms are proportional to the Yukawa couplings so that 
mixing effects can be neglected for the first two generations of sfermions. 

The lightest sneutrino can thus become the lightest SUSY particle (LSP) and a viable 
thermal DM candidate. Because of the large sneutrino mixing, the mainly RH sneutrino is 
no longer sterile, its couplings to SM gauge and Higgs bosons are driven by the mixing with 
its LH partner. Sufficient mixing provides efficient annihilation so that one can obtain a 
value for the relic density of Qh 2 ~ 0.11 as extracted from cosmological observations [OHH]. 

Direct detection (DD) experiments pose severe constraints on Dirac or complex scalar, 
i.e. not self-conjugated, DM particles because the spin-independent elastic scattering cross- 
section receives an important contribution from Z exchange, which typically exceeds ex- 
perimental bounds. In the mixed sneutrino model, this cross-section is suppressed by the 
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sneutrino mixing angle. Therefore, on the one hand a viable sneutrino DM candidate re- 
quires enough mixing to provide sufficient pair-annihilation, on the other hand the mixing 
should not be too large in order not to exceed the DD limits. Here we will explore the 
parameter space of the model where these conditions are satisfied. 

It is intriguing that a mixed sneutrino also opens the possibility for a supersymmetric DM 
candidate below 10 GeV. Light DM candidates have received a lot of attention recently [§]- 
[T9] because of results of DD experiments that show hints of events compatible with light 
DM. This includes the modulation signal from DAMA [20] as well as recent results from 
CoGeNT [21] and CDMS [22]. The best fit values for the mass and the cross-section do not 
overlap when these results are interpreted as a spin-independent contribution, nevertheless 
analyses taking into account the signal's dependence on the DM velocity distribution have 
shown [T2K231I21] that the observed events can be compatible with the null results obtained 
by other experiments such as Xenon [25j l26]. In addition, CRESST-II has very recently 
reported 32 events with an expected background of 8.7 ± 1.4, compatible with a DM mass 
of 15 GeV or below and a spin-independent cross-section of a few times 10 -5 pb [27]. 

Whether or not these events are confirmed, the possibility of light DM with large elastic 
scattering cross-sections remains interesting and particularly challenging to probe because 
experiments suffer from a severe loss of sensitivity at low masses. Furthermore, the mixed 
sneutrino revives the possibility of light DM in the MSSM: although it is possible to find a 
light neutralino that has the right properties to satisfy DM constraints in the MSSM with 
non-universal gaugino masses, this scenario is incompatible with additional constraints on 
the model [2"B~||2T?] . Indeed efficient annihilation requires additional light particles, for example 
a second doublet Higgs with mass in the 100 GeV range, which is strongly constrained by 
B-physics processes, in particular B s — > fi + n~ [30]. A light neutralino is still possible in 
singlet extensions of the MSSM where new Higgs singlets provide additional possibilities for 
efficient annihilation of the lightest neutralino [T8"|I3"T]. 

The phenomenology of the mixed-sneutrino model that we examine here was first investi- 
gated in [5]. Indirect detection signatures were discussed in [32], and LHC signatures in [33] PI 
We extend on these analyses in several ways. First of all, in contrast to the above mentioned 
studies, we here concentrate on light DM with mass of about 10 GeV and below. Second, we 
explore the parameter space of the model that gives a consistent sneutrino DM candidate 
using up-to-date constraints on elastic scattering cross-sections, examining also the effects 
of uncertainties in, e.g., the DM velocity distribution. In our scans, we take into account 
radiative corrections to the SUSY and Higgs spectrum; in particular we include the 1-loop 
corrections to the sneutrino mass originating from the Higgs contribution, and those to the 
light Higgs mass originating from the large A„ term. Moreover, we consider both the case 
of one and of three sneutrino flavours, assuming complete degeneracy in the three flavour 
case. For the allowed scenarios, we explore the consequences for DD as well as for indirect 
detection in photons, antiparticles and neutrinos. Finally, we explore the consequences for 
searches at the LHC and ILC. 

We characterize the scenarios that satisfy the DD constraints, including those that are 
within the region favoured by CoGeNT (and maybe also CRESST). The allowed scenarios 
have specific characteristics which include, e.g., dominantly invisible Higgs decays. Besides, 
if the charged sleptons are heavier than the X2 an d Xii as is the case over most of the valid 
parameter space, this implies dominantly invisible decays of neutralinos (x?2 ~~ v ^\) an d 

1 Many more studies of sneutrino DM have been performed in other models, like models with extra singlets 
or models with Majorana neutrino masses, see e.g., [TTll5H - fli?] . 
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single- lepton decays of charginos (xf — > i^vi). 

The paper is organized as follows. Section [2] describes the framework of our analysis, 
giving details on the model, the mass spectrum, and radiative corrections. Section |3] then 
discusses collider constraints, sneutrino annihilation and direct detection. The relic density 
and DD predictions for the one light sneutrino case are analyzed in Section HJ including a 
discussion of astrophysical uncertainties. Results for three degenerate light sneutrinos are 
presented in Section [5j Signatures in indirect detection are discussed in Section [6] and collider 
signatures in Section [71 A summary and conclusions are given in Section [8j The Appendix 
contains Feynman rules for the relevant sneutrino interactions. 

All numerical results have been obtained with micrOMEGAs [50|[5Tj . linked to an appro- 
priately modified version of SuSpect [52J. 



2 Framework 

2.1 Mixed sneutrinos 

The framework for our study is the model of [5] with only Dirac masses for sneutrinos. In 
this case, the usual MSSM soft-breaking terms are extended by 

A£ soft = m%, \Ni\ 2 + ApliNiHu + h.c. , (1) 

where m 2 ^ and A„ are weak-scale soft terms, which we assume to be flavour-diagonal. Note 
that the lepton-number violating bilinear term, which appears in case of Majorana neutrino 
masses, is absent. Neglecting the tiny Dirac masses, the 2x2 sneutrino mass matrix for one 
generation is given by 

M 2 = , m| + im|cos2/3 -^A vsm^ ^ 



■^Apvsin/3 



m\ 



Here m| is the SU(2) slepton soft term, v 2 = v\ + v\ = (246 GeV) 2 with u 1:2 the Higgs 
vaccuum expectation values, and tan/3 = Vijv\. The main feature of this model is that the 
m~, m 2 ~ and Ay are all of the order of the weak scale, and A v does not suffer any suppression 
from Yukawa couplings. In the following, we will always assume m^ < mi so that the lighter 
mass eigenstate, Ui, is mostly a. vr. It is in fact quite natural to obtain this relation when 
embedding the model in a GUT scale model, because the renormalization group running of 
mi is governed by M 2 , while for m^ the running at 1 loop is driven exclusively by the A v 
term, since Vr is a SM singlet. 

A large A v term in the sneutrino mass matrix will induce a significant mixing between 
the LH and RH states, 

v x = cos 9p v R - sin 9 U v L , (3) 
z> 2 = sin Op v R + cos 9 U v L , (4) 

where m^ < m„ 2 and the mixing angle 
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Notice that for fixed sin^, A u is proportional to m? 2 — m^. This means that A v is of the 
same order as other soft terms in the sneutrino sector. 

A large value of A v can induce a large splitting between the two mass eigenstates even 
if m| and m 2 ^ are of the same order, leading to scenarios where m Vl <C m U2 ,mj L . In this 
way, V\ can naturally be driven much below the neutralino masses. The model can easily be 
generalized to three generations. When doing so we will neglect for simplicity any flavour 
mixing in the sneutrino sector. 

The couplings of the mostly sterile v\ are those of the LH sneutrino suppressed by a 
factor sin 2 9 U due to mixing. In addition, there is a new direct coupling between the Higgs 
bosons and the LH and RH sneutrino components proportional to A v . The couplings to the 
Z and light Higgs boson will play a crucial role both for annihilation processes and for the 
elastic scattering cross-section. The coupling of sneutrinos to neutrinos and neutralinos is 
dominated by the wino component of the neutralinos and will be important for the annihi- 
lation of sneutrinos into neutrinos. The relevant Feynman rules are given explicitly in the 
Appendix. 



2.2 Particle spectrum 

We assume a model with soft terms defined at the weak scale and unification of gaugino 
masses at the GUT scale. The latter leads to M 2 ~ 2M\ ~ M 3 /3 at the weak scale. For the 
sneutrino sector, we take the masses and mixing angle, m Vl , m V2 and 6„, as input parameters, 
and compute the m^, and Ay. This also fixes the corresponding LH charged slepton 
mass term; for the RH one we assume m R = m^. Note that this choice has no effect as 
concerns DM properties, but can have implications for collider searches as will be discussed 
in Section [TJ 

In the one-generation case, we assume that only the tau-sneutrino is light and all others 
are heavy, with soft masses of 1 TeV. In the three-generation case, on the other hand, 
complete degeneracy between the slepton generations is assumed. 

For the squark sector, we assume a common soft mass nig = 1 TeV and take A t = — 1 TeV 
in order to avoid the constraint on the light Higgs mass. Other trilinear couplings for charged 
sparticles are neglected. The higgsino mass, pseudoscalar Higgs mass and tan (3 are also input 
parameters, we fix them to fi = 800 GeV, Ma = 1 TeV and tan /3 = 10. 

We use a modified version of SuSpect [52] for the spectrum calculation. The original 
SuSpect includes the 1-loop radiative corrections to neutralino, chargino and squark masses; 
corrections to Higgs masses are implemented at the two-loop level. We have extended it to 
include RH sneutrinos, and implemented 1-loop radiative corrections to sneutrino masses as 
well as those to the light Higgs mass induced by the A v term. 



2.2.1 Radiative corrections to sneutrino masses 

Let Ml(Q) be the running mass matrix Eq. ([2]) at the renormalization point Q. Then, the 
pole mass matrix is given by 

Ml = Ml(Q)+\ uLLW) -lrW) \ 6 

where AM~ LL (Q), AM~ LR (Q), AM? RL (Q) and AM~ RR (Q) are radiative corrections to the 
LL, LR, RL and RR components, respectively. For large Ay terms, contributions to AM|(Q)'s 
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Figure 1: Tree- level and 1-loop soft terms 
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are dominated by the one-loop slepton-Higgs diagrams. The main contribution to the light 
sneutrino mass shift results from AM? RR (Q) , which can be aproximated as 

AMf M (Q)| app = ^(log^-lJ • (7) 

Figure [TJ shows the the sneutrino soft terms mj, and that give a loop-corrected 
light sneutrino mass of 5 GeV and mixing angle sin 9o = 0.3 as a function of m„ 2 . The 
momentum and renormalization scales are set at p = Q = (m£; 1 m£; 2 ) 1 / 2 |i_i 00 p. The figure also 
compares the complete 1-loop result (dashed green lines) to the approximation m 2 ^ | treejapp = 
^i>i 1 1— loop — & m drr{Q)Upp (dotted blue lines). 

2.2.2 Higgs mass corrections 

After minimizing the Higgs potential, the mass of the lightest neutral Higgs reads 

m h = m l sin 2 (a + (3) + m\ cos 2 (a — (3) 

+ v 2 [(AAis 2 c 2 + AA 2 c 2 ^ - (AA 3 + A\4)c a s a cpSp 

+AA 5 (c 2 cJ + s 2 ^) - 2(AA 6 s aC/3 - AA 7 c aS/3 ) cos(a + (3)] (8) 

with s a = sin a, c a = cos a, etc.. The AAj include the radiative corrections to the quartic 
couplings. Loop diagrams involving sneutrinos can induce corrections to the quartic cou- 
plings through the presence of the weak scale A„ term. If we neglect the Yukawa couplings 
of the sleptons as well as the trilinear terms for the charged sleptons, then only X 2 receives 
a correction of 

1 IA,I 4 (™% +1^- m~ \ 
— 2 2,2 [-T (9) 

i=l v V2 v\' \ i>2 v\ v\ / 
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Figure 2: Contours of the radiatively-corrected light Higgs mass (from top to bottom: 111, 
114, 115 and 115.2 GeV) in the sin# 5 versus plane, for = 5 GeV. The dotted red 
lines are for the case of one light sneutrino, while the full black lines are for three light 
sneutrino generations. 

with the sum running over the sneutrino flavours, Nf = 3. Note that AA^ is negative, thus 
resulting in a decrease of the light Higgs mass Eq. [HJ We have checked that the effective 
potential technique gives the same result as PBMZ [53] when p 2 = 0. Setting p 2 = m\, we 
obtain the Higgs pole mass. We have incorporated the corrections to the light Higgs mass 
due to the sneutrinos in SuSpect. When A„ is large (which means large mixing and a large 
m„ 2 ) the radiative corrections from the sneutrino sector can drive the light Higgs mass below 
the LEP limit. This is illustrated in Fig. El where we show contours of constant rrih in the 
sinOp versus m^ 2 plane. The dotted red lines are for the case of one light sneutrino, while 
the full black lines are for three degenerate light sneutrinos; in either case, mo 1 = 5 GeV. 
Notice that the contour of = 111 GeV for three contributing sneutrino flavours almost 
falls together with the = 114 GeV contour of the one-sneutrino case. 

3 Constraints on the model 
3.1 Collider constraints 

A light sneutrino with nip < will contribute to the invisible width of the Z boson, thus 

putting a constraint on the sneutrino mixing: 

N f ■ An ( ,o x2\ 3 / 2 

^gr„^(l-(^)) <2MeV (10) 

where T u = 166 MeV is the partial width into one neutrino flavour. For one light sneutrino 
with = 5 (20) GeV, this leads only to a mild constraint on the mixing angle of sin 9$ < 
0.39 (0.43). In the case of three degenerate sneutrinos, this constraint becomes stricter, 
sin0p < 0.296 (0.33). 
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We also impose the limits from SUSY [51] and Higgs [55J searches at LEP2. Accounting 
for a theoretical uncertainty in the light Higgs mass of about 3 GeV, we require rrih > 
111 GeV. For a large value of the sneutrino mixing this implies an upper bound on m U2 , see 
Fig. |2j The radiative processes where a photon is emitted in addition to a pair of invisible 
supersymmetric particles will contribute to the process e + e~ — > 7 + invisible, which has 
been searched for by the LEP2 experiments. Here invisible particles include not only the v\ 
but also Xi or even when they are the NLSP and NNLSP respectively as they decay in 
v\v. As the LEP2 limit we take e + e~ ->■ 7 + invisible < 0.15 pb at y/s = 189 - 209 GeV 
for p\ > 0.02 y/s and ^ eam7 > 14 deg [56J. We have computed the full 3-body cross-section 
for the single photon production using calcHEP [57] and found that it rarely exceeds tens 
of fb, which means it does not constrain our model. The reason for this is that for a 
large cross-section it is necessary to have a light particle exchanged in the t-channel. For 
neutralino production this means a light selectron, which is only possible in our 3 generation 
model. Likewise, VeVel production, which is enhanced by t-channel chargino exchange, is 
contributing only in the 3 generation case. 

In [5B1G3] ^ was argued that the search for one-jet events with large missing transverse 
energy, so-called monojets, could provide a stronger limit on light DM than current DD 
experiments. Monojet searches at the Tevatron look for events with leading jet pr > 80 GeV 
and missing Et > 80 GeV, while 2nd jet px < 30 GeV and more jets are vetoed. An 
analysis [60] by the CDF collaboration of lfb -1 of data gave 8449 events, with an expected 
background of 8663 ± 332. Using calcHEP, we have computed the cross-section of pp — > 
Vi&i+g at Tevatron energies and found that after cuts it is typically of the order of 0.1 — 1 fb. 
Indeed for the parameter points that pass all other (including DD) constraints, see the scan 
of the following section, we find cross-sections of at most 1.5 fb. Thus the monojet search 
does not provide any additional constraint on the model. 

In what follows, when we discuss DM allowed scenarios, it is implicitly understood that 
collider constraints are satisfied. 

3.2 Relic abundance of sneutrino 

For computing the sneutrino relic abundance, we assume the standard freeze-out picture. 
We do not consider non-thermal sneutrino production. This is justified because the mixed 
sneutrino has electroweak interactions. We have implemented the mixed sneutrino model 
in micrOMEGAs 2.4, which allows for a fully automatic computation of the annihilation and 
DD processes. Note that in the computation of the relic abundance we have not included 
the extra degrees of freedom corresponding to the RH neutrino. As these particles decouple 
early, this will only induce a correction at the few percent level on the effective degrees of 
freedom, which is negligible for our purpose. 

The main annihilation channels for a light sneutrino are i) V\V\ — >■ vv (P*v* — > vv) 
through neutralino t-channel exchange, ii) V\v\ — > bb through exchange of a light Higgs in 
the s-channel, and Hi) V\v\ — > ff through Z exchange. The annihilation into neutrino pairs 
proceeds mainly through the wino component of the neutralino and is proportional to sin 9 U ; 
it is largest for light winos. The Z exchange is also proportional to sin 4 #,>. The light Higgs 
exchange, on the other hand, is proportional to (A u sin^) 2 . Note in particular that for a 
fixed value of the sneutrino mixing angle, the Higgs contribution will increase with as 
Ay also increases. 

The behaviour of Qh 2 as a function of the sneutrino mass and mixing angle is displayed 
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Figure 3: In black (full lines), contours of Qh 2 = 0.1 in the sin 6^ versus mo 1 plane for 
one light sneutrino, with = (200, 500) GeV and M 2 = (105, 200) GeV. The remaining 
parameters are fixed as explained in Sec. 12.21 The dashed red and dash-dotted blue lines 
show contours of constant ajyy = 10~ 4 pb and 10~ 5 pb, respectively. The dotted green line 
shows the limit from the Z width. 



in Fig. |3]for the case of one light sneutrino. A larger mixing is required for light masses. This 
is related to the fact that Qh 2 is inversely proportionnal to the number of degrees of freedom 
(g ff)- At the temperature where the QCD phase transition occurs, around Tqcd ~ 300MeV, 
the number of degrees of freedom starts to drop and Qh 2 increases. This is relevant for DM 
masses below ca. 6 GeV, where the freeze-out temperature Tf « ??7,dm/20 is of the order of 
Tqcd- Furthermore note that the uncertainty in the change of g e s around Tqcd will induce 
some uncertainty in the computation of Qh 2 . This is particularly important for < 2 GeV 
because in this case Tf m 100 — 150 MeV, precisely where there is a sharp drop and a large 
uncertainty in g e ^. We have not considered these corrections to the relic abundance as only 
a few scenarios fall in this category. 

The dependence of Qh 2 on the gaugino mass, M 2 is also displayed in Fig. |3j A lower 
mass requires a smaller mixing, this is because the self-annihilation channel into neutrinos 
is increased in this case. The relative contributions of the various annihilation channels are 
shown in Fig. H]for the two cases of ay large mixing, sin#,> = 0.35, but relatively small A„ 
(left panel) and b) smaller mixing, sin#;> = 0.22, but large A„ (right panel). In the large 
Ac, case, the contribution from the Higgs exchange enhances the bb channel if kinematically 
accesssible. 

3.3 Direct detection 

The spin-independent (SI) scattering of v\ on nucleons occurs through Z or Higgs exchange. 
The Z exchange is again suppressed by the sneutrino mixing angle. The scattering cross- 
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Figure 4: Relative contributions of different annihilation channels as a function of the v\ 
mass, on the left for = 200 GeV and sin^ = 0.35, on the right for = 500 GeV and 
sin 0c = 0.22, cf. Fig. EJ In both plots, M 2 = 200 GeV. 



section on a nucleus due to Z exchange is given by 
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4In = 1^/4 {(A - Z) - (1 - 4 sin 2 6 w )Zy sin 4 9 V , (11) 

where \i x is the sneutrino-nucleus reduced mass, while A is the atomic weight and Z the 
number of neutrons of the nucleus. One peculiarity of the Z-exchange contribution is that 
the proton cross-section is much smaller than the neutron one, with the ratio of amplitudes 
fp/ fn = (1 — 4 sin 2 9w)- The Higgs contribution on the other hand, which becomes dominant 
for large values of Ap, is roughly the same for protons and neutrons, 

^ = i ifet ( {A ~ Z)S £ 9h ™ f > n ' ( 12 ) 

where gh qq = e/ {2M w sw)x q (x u = — cos a/ sm(3,Xd = sin a/ cos/3) is the Higgs coupling to 
quarks after the quark mass has been factored out, and gho^ is the coupling to the LSP as 
given in Appendix A. 

The total SI cross-section is obtained after averaging over the v\N and v^N cross sections, 
where we assume equal numbers of sneutrinos and anti-sneutrinos. Here note that the 
interference between the Z and Higgs exchange diagrams has opposite sign for i>iN and 
v^N, leading to an asymmetry in sneutrinos and anti-sneutrinos scattering if both Z and 
Higgs exchange are important. We will come back to this in Section 16.11 when we discuss 
signals from sneutrinos captured in the Sun. 

A comment is in order concerning theoretical uncertainties. The computation of the elas- 
tic scattering cross-section for the Higgs-exchange diagram depends on the quark coefficient 
in the nucleons, which can be determined from the pion-nucleon sigma term, a^N and from 
the SU(3) symmetry breaking effect, uq = 35 ±5 MeV |61j . By default we take do = 35 MeV 
and a n N = 45 MeV [6 1 J . This leads to 

f d = 0.026, /£ = 0.020, = 0.13, 

P d = 0.036, # = 0.014, # = 0.13. (13) 

More recent estimates of the pion-nucleon sigma term typically indicate larger values of 
a vN = 55 — 73 MeV [52]- The very recent lattice results also tend towards a larger value 
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for the strange quark content of the nucleon, although uncertainties are still large 
The overall theoretical uncertainty that arises from the uncertainty in the scalar coefficients 
is relevant only for cases where the Higgs-exchange contribution dominates, since for the Z 
contribution the vector coefficients are simply determined by the valence quark content in 
the nucleon. Since the values in Eq. ( lT3"j) are rather on the low side, the DD cross-section 
shows a larger upward than downward variation when a n N and o"o are varied within their 
allowed ranges. For example, for (ovjv, o"o) = (70, 30) MeV the cross-section can increase by 
up to a factor 3.5 when Higgs exchange dominates, while for (55, 35) MeV the increase is 
at most a factor 1.7. The choice (45, 40) MeV on the other hand implies a decrease in the 
cross-section that can reach 30%. 

The limits on a^ 1 from DD experiments are extracted from the observed limit on the LSP- 
nucleus scattering cross-section assuming that amplitudes for protons (f p ) and neutrons (f n ) 
are equal. In our model this is not the case when Z exchange dominates. Therefore, we 
compute instead the normalized cross-section on a point-like nucleus, 



41* = ^ JP 1 v " ~ Jjnj (14) 



^ 2 x (Zf p + (A-Z)f, 



7T 



A 2 



where the average over v\ and v\ is assumed implicitly. This cross-section can be directly 
compared with the limits on a? 1 given by the experiments. For Xenon A = 131, Z = 54 
while for Germanium A = 76, Z = 32. 

The SI scattering cross-section mainly depends on the three parameters of the sneutrino 
sector as long as is large enough so that heavy Higgs exchange can be neglected. Con- 
tours of constant cr-Jjy in the sin 9 U versus m Ul plane are displayed in Fig. [3] for two different 
values of m V2 . For m/> 2 = 200 GeV, in the region > 6 GeV the A„ parameter is so 
small that scattering proceeds mainly via the Z-boson exchange; the detection rate is almost 
independent of . For m„ 2 = 500 GeV, Ay is so large that the nucleon scattering has an 
important component from the Higgs-boson exchange; for this channel the detection rate is 
proportional to Therefore the detection rate gets smaller for larger v\ mass. This is a 

characteristic of scalar dark matter particles. Similarly for m vx < 6 GeV, the nucleon scat- 
tering proceeds mainly through Higgs-boson exchange implying an increase for lower v\ mass. 

Before we proceed, another comment is in order concerning the tan /3 dependence. So far 
we have only considered tan/3 = 10. The main effect of increasing tan /3 is a lower f\ mass 
for the same choice of m Vl , , sin Qy and fi. In our approach, this can be translated into 
an upper limit on tan /3 as a function of m V2 , or vice- versa a lower limit on m U2 as a function 
of tan /?. In fact, for m V2 = 200 GeV like in Figs. [3]andHl one can only go up to tan /3 ~ 20; 
for higher values one first violates the fi mass bound and then gets tachyonic staus. For 
high values of, e.g., tan /3 = 50, one needs m V2 > 300 GeV. Another effect is a slightly higher 
h° mass for higher tan/3. 

The influence on the relic density of the z>i is small, the main effect being a slight suppres- 
sion of the V\0\ — > t + t~ channel through the interference of the Higgs and chargino exchange 
diagrams. To give a concrete example, for m Ul = 10 GeV, m U2 = 500 GeV, sin^ = 0.22, 
M 2 = 200 GeV and tan (3 = 10 we have flh 2 = 0.12 with 6% contribution from v x v\ — > t + t~, 
cf. FiglU For the same parameters but tan /3 = 50, we get Qh 2 = 0.131 with 4% contribution 
from v x v\ t + t~. Lowering M 2 to 105 GeV leads to Vth 2 = 0.092 (0.1) for tan/3 = 10 (50). 

The SI scattering cross-section is affected by tan /3 through the Higgs mass and the bottom 
Yukawa coupling, but this effect is also small. In the above example with M 2 = 200 GeV, 
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af ip = 8.63 x 1(T 5 pb and af rn = 1.17 x 10^ 4 pb change to af ip = 8.04 x 10~ 5 pb and af in = 
1.11 x 10~ 4 pb when increasing tan/3 from 10 to 50. Here note that rrih = 114.2 (115.2) GeV 
and Aj> = 310.7 (309.2) GeV at tan/3 = 10 (50). 

We conclude that the choice of tan (3 has little influence on the results. However, lower 
tan/3 has more parameter space in the sense that one can go to lower z> 2 masses. Therefore 
in the following we will focus on tan /5 = 10. 

4 Results for one light sneutrino flavour 

For the light sneutrino to be a viable dark matter candidate, in addition to the collider 
constraints, we require that both the 3a upper bound from WMAP5 [B], flh 2 < 0.1285, and 
the DD limits be satisfied. In fact, DD here provides the more stringent constraint. We first 
discuss the case of one light sneutrino, concentrating on the third generation. 

For the numerical analysis, we perform a random scan of the parameters of the sneutrino 
and gaugino sector. To search efficiently the region with a light sneutrino, we use mp 1 2 and 
sin#p as input, from which we compute mj-, and = A„ 3 . The parameter ranges we 
scan over are 



All other parameters have little influence on the light sneutrino scenario and are therefore 
fixed as specified in Section 12.21 In particular, the stau masses are determined through 
m R3 = m L 3 an d A T = 0, while all other soft masses are set to 1 TeV. Moreover, having 
checked that it does not influence the results, we set fi = 800 GeV. The scan is carried out 
in two steps with 50000 points each, the first varying M 2 from 100 to 150 GeV (light wino 
case) and the second varying M 2 from 150 to 500 GeV (heavy wino case). Extending the 
upper value for the range of M 2 has no influence on our results. 

The scan results are shown in Fig. [5] in the sin 0y versus and Qh 2 versus planes. 
We find scenarios that satisfy the WMAP and DD limits for sneutrinos as light as 1 GeV. 
Such light 2/1's below ca. 4 GeV require a large mixing, sin6^ > 0.25, and annihilate predom- 
inantly into neutrinos. A more modest mixing of sin 9$ > 0.12 is needed for V\ masses above 
the 6-threshold, where annihilation into bb through Z or h can contribute significantly. 

Figure [6J shows the expected SI cross-sections for Xe as a function of the i>\ mass, to- 
gether with the limits from CoGeNT, XenonlO and XenonlOO, which give the most stringent 
constraints on light DM [2TJ[25J[25|. The cross-section is appropriately re-scaled if the V\ is 
only a part of the DM, i.e. £ = f2/i 2 /0.11 for points with Qh 2 < 0.0913 and £ = 1 otherwise. 

When the sneutrino annihilation is dominated by Z or Higgs exchange, the prediction 
for the SI scattering cross- sect ion, which is also dominated by the Z or Higgs exchange 
diagrams, is directly related to the annihilation cross-section. In this case, after re-scaling 
for the lower dark matter density if Qh 2 < 0.0913, crf^ varies only within a factor of 2 (not 
taking into account uncertainties due to the quark coefficients Eq. (|13p ). On the other hand, 
when the sneutrino annihilates dominantly into neutrino pairs, there is no such correlation 
between Qh 2 and <jf* N , and it is possible to suppress the direct detection cross-section by 
more than one order of magnitude. These scenarios correspond to the scatter points with 
low cross-sections in Fig. 



1 GeV < m 0l < 15 GeV, 
< sin 5 < 0.5 , 



100 GeV < m V2 < 1000 GeV , 

100 GeV < M 2 = 2Mi < 500 GeV . 



(15) 
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Figure 5: Scan results for the one light sneutrino case. On the left, a scatter plot of sin 6$ 
versus mo 1 - Here the red (blue) dots show models that have a relic density within (below) 
the 3cr WMAP range. The yellow triangles refer to models that satisfy the WMAP upper 
bound as well as the DD limits from CoGeNT and XenonlO. The plot on the right shows 
Vth 2 versus with the light blue points passing the DD constraints. 

The points that successfully pass all constraints have sneutrino masses of about 1 to 
8 GeV and include many scenarios that have cross-sections within the range favoured by 
CoGeNT. In particular the bulk of our scenarios with masses around 7-8 GeV are in the 
region allowed by CoGeNT and Xenon. Recall also that such light DM with a cross-section of 
a few times 10 -5 pb could cause the events observed by CRESST. For mp 1 ~ 5 GeV, many of 
the WMAP-allowed scenarios lie above the CoGeNT limit. A value of M 2 ~ 100 - 120 GeV 
is required in this case to achieve efficient enough annihilation into neutrinos. For lower 
masses, the DD limit becomes much weaker, leaving the relic density and Z width as the 
main constraints. Moreover, we recall that, while some fine-tuning is needed to achieve a 
light there are no large hierarchies among the soft terms (cf. Section |2"TT|) . In particular 
the ratio between Ay and M 2 roughly ranges from 1/5 to 5 (8) for < 6 (8) GeV. 

It is also intriguing that in Fig. [6] there is a lower limit on the DD cross-section of about 
10~ 5 pb, which is almost flat for mo 1 > 4 GeV. It arises because the lowest oJFjy is obtained 
when the dominant contribution to the relic abundance is annihilation into neutrinos. In this 
case, scattering off nucleons proceeds dominantly through the Z-boson exchange, which is 
almost independent of m Sl . Note that this lower limit is free of the theoretical uncertainties 
on the quark coefficients in the nucleon. The mild mass dependence of the lower limit on 
erf ^ N for Tdvy < 4 GeV arises because the Higgs exchange which is dominant in this region, is 
porportional to l/m~ , Eq. (TTTT) . Overall, this lower limit means that mixed sneutrino DM 
with a mass of a few GeV will be either discovered or excluded if the DD experiments can 
cover SI cross-sections down to 10~ 5 pb. 

The mass of the z> 2 and consequently of the fi 2 is also constrained. When M 2 is small, 
around ca. 100 GeV, and annihilation into neutrinos dominates, the angle sm9a > 0.25 to 
provide sufficient annihilation but the mass of z> 2 is not constrained. On the other hand when 
Higgs exchange is dominant, the annihilation is proportionnal to A\ sin Q\ m m% sin 9~, which 
gives a lower bound on . Consequently a lower bound on m fl > 250 GeV is obtained for 
large M 2 . This is illustrated in Fig. which shows the points of Fig. [5] in the m-± versus m f 

Al 
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Figure 6: Scatter plot of re-scaled a~^ Xe as a function of for the case of one light sneutrino. 
The red (blue) dots show models that are within (below) the WMAP range. Also shown 
are the limits from CoGeNT (green), XenonlO (dashed), and XenonlOO limit (crosses). The 
CoGeNT and XenonlO limits are for vq = 220 km/s, v esc = 600 km/s and p = 0.3 GeV/cm 3 . 
The XenonlOO band is obtained by varying vq = 180-260 km/s, v esc = 500-600 km/s, 
p = 0.2-0.6 GeV/cm 3 and an L eff that decreases at low energies, see Sec. 14.11 for details. 

plane. We found a few scenarios with < rn^±. For none of these points X2 or Xi have 
a significant branching fraction into slepton/lepton pairs. This has important consequences 
for collider searches, as we will discuss in Section [7J 

Let us briefly come back to the tan /? dependence. We have argued above that tan /3 = 10 
has a larger parameter space than higher values, but otherwise results are very similar. This 
is illustrated in Fig. [HI which shows scan results as in Fig. [5] but for tan/3 = 50. As can be 
seen, the only change is that points with low DD cross sections, corresponding to a low , 
are removed. Indeed the scan points for tan (3 = 50 that pass all constraints are just a subset 
of those for tan /3 = 10. In the following we will therefore use the tan /3 = 10 scan points. 

4.1 Uncertainties in the DD limits 

As the DD limit plays a crucial role in constraining the parameter space, it is important to 
consider the uncertainties involved in extracting the limit on a~£ from experiments. One 
source of uncertainty is the assumed DM velocity distribution. In particular, there is a 
minimum velocity necessary to pass the energy threshold of the detector, v 2 min = ^0^. Thus 
the direct detection limit for low masses is particularly sensitive to the upper part of the 
velocity distribution. To illustrate this effect, we assume an isothermal velocity distribution 
and allow for a 2a variation of vq = 220 ± 20km/s. We also vary the velocity above which 
DM will escape from the Galaxy in the range 500 < v esc < 600km/s [65] H The local dark 
matter density is another factor that possesses a large uncertainty, it is generally assumed to 

2 Several velocity distributions were considered in [661167] . Varying the input parameters of the isothermal 
distribution as we do, reproduces most of the span of variation of the limits found in those studies. 
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Figure 7: Same as the lhs plot of Fig. |5] but in the m~± versus plane. The red (blue) 
dots denote points that have a relic density within (below) the 3a WMAP range but too 
large a DD cross- sect ion. The yellow triangles denote points that satisfy the WMAP upper 
bound as well as the DD limits from CoGeNT and XenonlO. 

be p = 0.3 GeV/cm 3 but can vary from 0.2 — 0.7. (For a review of astrophysical uncertainties 
on the velocity and density distribution see 

Finally, there is an additional uncertainty that is specific to Xenon, it comes from the 
scintillation factor L c g that allows to convert the measured electron energy into the nuclear 
recoil energy. In particular the threshold for nuclear recoil energy strongly depends on the 
scintillation factor at low energies. New measurements have shown that the scintillation 
efficiency can be significantly lower than the value used in Ref. [25] — and this especially at 
low energies — thus weakening the limit on [HSj- The XenonlO limit in Fig. uses a 
constant value for L cff while the band for XenonlOO uses a value of L eS that decreases 
at low energies [70] . 

5 Results for three degenerate light sneutrinos 

In the case of three light sneutrinos, the constraint from the Z width on the sneutrino mixing 
angle is stronger. Moreover, for three exactly degenerate sneutrinos, the relic density can in- 
crease by a factor up to 3 as compared to the one-generation case. This is because the effective 
annihilation cross-section including the coannihilation channels is (av) oc (^ df^ij)/ XI 9? 
where gi is the number of degrees of freedom and the cross-section for annihilation of two 
supersymmetric particles, Xii Xj SM particles. Thus when annihilation is dominated by 
the processes ViV* — » //, (erf) increases by a factor of 3 as compared with the one genera- 
tion case. This then leads to a tension between the relic density constraint which requires 
significant mixing and the Z width which strongly constrain this mixing angle. The relic 
density constraint is particularly severe at low masses where we had Qh 2 »s 0.11 in the one 
generation case. The smaller mixing can of course be partially compensated by Ay however 
this then typically violates the DD limits. Furthermore, for three light sneutrinos with large 
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Figure 8: Same as Fig. |5]but for tan/3 = 50. 




2 4 6 8 10 12 14 16 



v 1 mass (GeV) 

Figure 9: Scatter plot of crfl Xe versus mo 1 analogous to Fig. |6]but for three degenerate light 
sneutrinos. 

A&, the Higgs mass constraint becomes much more severe. 

Decreasing makes it easier to satisfy the relic density constraint. Indeed in this case 
annihilation is dominated by production of neutrinos and all channels ViVj — > ViVj contribute 
equally to the annihilation cross-section so that Vth 2 is only a factor 3/2 larger than for the 
one generation case. 

We perform a random scan analogous to the one in the previous section. The result of 
this scan is shown in Fig. [9j As expected, the parameter space for sneutrino DM is now much 
more restricted, and only few points pass all constraints. The allowed sneutrino masses now 
range from 4.5 to 8 GeV. For most scenarios the chargino mass lies just above the LEP limit. 
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The sneutrino mixing is around sin^ ps 0.27 — 0.29, while the slepton masses are restricted 
to mi L ~ mp 2 « 150 — 480 GeV. This is because larger masses induce too large corrections 
to nih, while smaller masses reduce the value of Ap and thus the Higgs-exchange contribution 
to the annihilation processes. It is also interesting to note that in the DD-allowed region 
there are no points with VLh 2 below the 3a WMAP range (and overall there are only very 
few points with Vth 2 < 0.0913). 

The exact mass splitting between sneutrinos of different flavours strongly influences the 
allowed parameter space. For example, a 1 GeV mass splitting is enough to suppress any 
coannihilation contributions to a negligible level. Such a mass splitting can be induced by 
small splittings in the soft terms, which are rather generic in supersymmetric models even if 
one starts out with universal soft terms at a high scale. Taking, for instance, = = 
m Z 3 > m iVi = m N 2 = m N 3 an d J us t 1% difference in the Ap terms, Ap eli = 0.99 A„ T , we find 
mass splittings of 2-14 GeV. In this case, the heavier v\ e ,^ decay into the v\ T LSP through 
3-body modes [71] with decay widths ranging from 10~ 12 to 10~ 8 GeV. The dominant mode 
is into neutrinos; visible decays, e.g., v\ e — > e =F r ± i> lr , have at most few percent branching 
ratio. In this case we are back to the case of one light sneutrino presented previously with 
only a more severe constraint on sin Op and Ap from the Z width and light Higgs mass. In 
what follows we will therefore concentrate on the case of one light sneutrino. 

6 Indirect detection 

6.1 Neutrinos from annihilation in the Sun 

Neutrinos originating from annihilation of DM captured in the Sun could provide a good 
signature for the mixed sneutrino model. In scenarios with light sneutrinos, a detector with 
a low threshold is needed. The Super-K detector has a threshold of 1.6 GeV, while large 
detectors like Antares or ICECUBE have threshold energies above 25 GeV, making them 
unsuitable for detecting the neutrino flux from light sneutrinos. 

The capture rate for DM particles in the core of the Sun depends on the DM-nucleus 
scattering cross-section, as well as on the DM velocity distribution and local density. After 
being captured, the DM annihilates into Standard Model particles, which further decay into 
neutrinos that can be observed at the Earth. The capture rate is approximated as [72"|I73] 

where m^, the mass of the nuclear species i, and mp 1 are given in GeV. v is the DM velocity 
dispersion, is the mass fraction of element i in the Sun, and <pi its distribution. Fi(m) is a 
form factor suppression and S a kinetic suppression factor. For these parameters, we use the 
values listed in Tables 8 and 9 of Ref. [72]. Finally, crf^ is the elastic scattering cross-section 
on point- like nucleus. 

In models where the DM is not self-conjugate, one can get different capture rates for 
particles and antiparticles. Furthermore both particle-particle (antiparticle-antiparticle) 
and particle-antiparticle annihilation channels exist, in our case v\0\ — > vv (Akx) anc ^ ^i^i ~~ ^ 
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XX (A xx ). The equations describing the evolution of the number of DM (anti-)particles, 



N X (N X ), are then 



*x = C x 

The annihilation rates can be approximated as 



2A XX N 2 X - A XX N X N X , 



A X *N X N X , 
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(16) 



(17) 



where V P 



5.8 x 10 30 cm 3 (m x /GeV) 3//2 is the effective volume of the core of the Sun 



'eff ~ - — v ,o x/ 

If the capture and annihilation rates are sufficiently large, equilibrium is reached and the 
annihilation rate is only determined by the capture rate. We assume this to be the case. To 
take into account the different capture and annihilation rates for particles and antiparticles, 
we define 
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Then, after equilibrium is reached, N x = N x = and we can solve Eq. [T6j The annihilation 
rate at present is determined by the capture rates as well as a: 
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(20) 



The total neutrino spectrum at the Earth, assuming self-annihilation channels are solely into 
neutrino pairs, is given by 



J-[r Br ^ 



4vrrf2 \ L ™* r ™ dE 




(21) 



where d = 1.5 x 10 8 km is the distance from the Sun to the Earth, Br vv is the branching 
fraction for annihilation into neutrino pairs Br^j the branching fraction into each parti- 
cle/antiparticle final state //. Nf and N vv {Nvp) are the neutrino spectra resulting from 
those annihilations. Here dN uu /dE is simply proportional to a delta function. The neutrino 
spectrum originating from different annihilation channels into SM particles and taking into 
account oscillations was computed in [75], we use the tables given there. Note that for the 
neutrino pair an average over the three flavours in the annihilation process is assumed. For 
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the one light sneutrino scenario this is not the case, however this is still a good approx- 
imation since almost perfect 3-generation mixing is expected for neutrinos below 10 GeV 
propagating in the Sun [75] . 

For light DM particles, the process of evaporation from the Sun can be important, this 
effect was estimated in [731171] . We modify Eq. [TB] accordingly and solve it iteratively. We 
find that for the range of annihilation cross-sections of our scenarios, the evaporation affects 
significantly DM particles lighter than 3 GeV and is irrelevant for heavier DM particles. 

Finally, to compare with the data, one must compute the muon flux for upward events [76] 

m H oo E v , 

% - 1 dE -M S is S < I iikf Lp -^ m - EAK , ,» (22) 



where the survival probability for a muon of energy E' and final energy E^ 




and the muon energy lost after propagating a distance z is 

= - H(i _ e -ft») (24) 

Here y = m^/^crap), r is the muon lifetime, p = 2.6 g/cm 3 the rock density and a = 
2 x 10 -3 GeVcm 2 /g, (3 = 3 x 10~ 6 cm 2 /g characterize the average energy loss of the muon 
traveling through rock or water. dPcc is the probability for a neutrino with energy E v to be 
converted into a muon of energy E^ over a distance dz through charged current interactions. 

The event rate takes into account the effective area of the detector. For Super-K, the 
detector is cylindrical with a radius R = 18.9 m and a height 36.2 m. The muons that are 
stopped in the first 7 m are not observed, this corresponds to the muon energy threshold 
of E = 1.6 GeV. Furthermore, only muons of energy larger than 7.7 GeV go through the 
detector. This means that for our allowed scenarios with DM masses below 8 GeV, almost all 
the muons will be stopped within the detector. We therefore compute the rate for stopped 
muons only. To do this we take into account the zenith angle 8 Z when computing the effective 
area of the detector and average over — 1 < cos 6* < [77]. 

In Fig. [TU], we display the total muon/antimuon rate as a function of the v\ mass for the 
scan results of Section 0] (only the points that pass all constraints are used). The Super- 
K limit as extrapolated from [78] is also displayed. Here we take v = 270 km/s for the 
DM velocity dispersion and p^ — £ x 0.3 GeV/cm 3 with the same re-scaling factor £ as in 
Section H] when the sneutrinos do not make up all the DM. Furthermore, we take the values 
of Eq. f[T3|) for the quark coefficients. The larger rates are found for sneutrinos lighter than 
nit, that annihilate dominantly into t + t~ or neutrino pairs, since these modes give a harder 
neutrino spectrum. Some of our scenarios are within 10% of the experimental limit in the 
mass range 3 — 4 GeV. In general, the expected rates are within one order of magnitude of the 
Super-K limit, except for < 3 GeV where the evaporation effect becomes very important 
so that the neutrino flux is strongly suppressed. Note that these predictions depend on the 
DM velocity and average density. In particular a lower v would lead to larger neutrino rates. 
At the same time a lower DM velocity would relax the DD constraints. 
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Figure 10: Predicted flux from stopped muons from the Sun in the Super-K detector for 
v = 270 km/s (left) and the corresponding ratio of antimuon to muon fluxes (right), for 
the allowed scan points of Section 21 cf. Figs. |5] and [6j The red (blue) points have a relic 
density within (below) the 3a WMAP range. The horizontal line in the lhs plot shows the 
extrapolated Super-K limit. 



As already mentioned, one characteristic of our mixed sneutrino DM scenario is that the 
scattering rates on nucleons can differ significantly for particles and antiparticles. In fact, 
the capture of v* is in general more efficient than that of v due to the destructive interference 
between the Z and Higgs exchanges in the latter case. As a result, the flux for antineutrinos 
and thus of antimuons is often much larger than that for muons, see the rhs plot in Fig. [T01 
Distinguishing muon from antimuon events would therefore provide an additional test of this 
model. 

6.2 Photons 

As mentioned previously, when > 5 GeV the dominant v-y pair annihilation channels 
are v T v T or bb, while for lighter DM the charged fermion channels are cc and t + t~, with 
both channels having similar rates. Below ca. 1.5 GeV, the annihilation is purely into 
neutrinos. The annihilation channels into charged fermions leave a signature in photons, 
antiprotons and positrons. Photons are particularly interesting as the Fermi-LAT satellite 
is currently taking data in this channel. (The signature in charged cosmic rays is discussed 
in the next subsection.) The Fermi-LAT collaboration has obtained its first limits on the 
flux of photons originating from DM annihilations. In particular limits on av were extracted 
from observations of dwarf galaxies [79]. The best one is obtained from Ursa Minor and 
corresponds to av > 7 x 10~ 26 cm 3 s _1 for todm = 10 GeV assuming an annihilation entirely 
into bb. 

To study the sensitivity to light sneutrino DM, we compute the annihilation cross-section 
into the // channel, avjf, for the parameter points of Section |4] that pass all constraints. 
To this aim we assume a NFW dark matter profile [80]. The results are shown in Fig. [Til for 
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Figure 11: Annihilation cross-sections crv^j for / = b, r as a function of m^, for the allowed 
scenarios of Section HI 

f = b and r. The cc channel gives similar results to t + t~. 

When bb is the dominant charged particle final state, av h i ranges from 10 _27 cm 3 s _1 to 
10~ 25 cm 3 s~ x . We therefore expect that some of our scenarios could be probed by Fermi-LAT 
once the analysis is extended to lower masses. When t + t~ is the dominant channel, av Tf 
can reach up to lCT 26 cm 3 s -1 . However, this channel gives a harder photon spectrum, so 
that the sensitivity on the photon flux is expected to be almost an order of magnitude better 
than for the bb channel [79] . The Fermi-LAT satellite could therefore also probe some of the 
scenarios with i>\ masses below ca. 5 GeV. Note that for < 2 GeV, the annihilation is 
dominated by neutrinos, leaving no signature in the photon channel. These scenarios are 
also the ones that because of evaporation have very suppressed rates in neutrino telescopes. 

6.3 Charged particles: positrons and antiprotons 

The annihilation into charged particles will also leave a signature for antiprotons and positrons 
The flux of antiprotons has been measured by PAMELA [81] together with the ratio of an- 
tiprotons to protons. It is well described by the expected background flux from standard 
astrophysical processes, the spallation of cosmic ray proton and helium nuclei over the in- 
terstellar medium. There are, however, large uncertainties in the theoretical predictions 
of the secondary fluxes, due to uncertainties in the p production cross-sections, as well as 
in the parameters of the propagation models (82]. Each of these uncertainties have been 
estimated to be around 25% for the diffusion model with convection and re-acceleration 
described in [S2]. Therefore there is still room for an additional contribution from DM 
annihilation, especially in the lower part of the spectrum where the experimental uncer- 
tainties are the largest. For an energy of Ep = 0.56 GeV, for instance, which roughly 
corresponds to the peak of the spectrum from a 5 GeV DM particle, the measured flux is 
$p = 15.3j£5 ± 0.9 x 10~ 7 (GeV cm 2 ssr) 1 . The left panel in Fig. [121 shows antiproton 
spectra originating from V\ annihilation for = 1.4, 3.9 and 7.5 GeV from left to right 
for either line-style. For each of these masses we have selected two scenarios corresponding 
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Figure 12: On the left, the antiproton spectrum for 6 representative allowed points as ex- 
plained in the text; the antiproton flux measured by PAMELA [HI] is also displayed (crosses). 
On the right, the differential flux for V\ annihilation into antiprotons at E = O.lm^ includ- 
ing propagation for two sets of propagation parameters, MIN (light blue points) and MED 
(dark blue points), see Table [TJ 

to near maximal (full lines) and near minimal (dashed lines) flux. The selected scenarios 
are within 10% of the absolute maximum/minimum fluxes. As can be seen, the predictions 
can vary by around two orders of magnitude for a given mass, depending on the details of 
the parameter point. Indeed, while the antiproton fluxes from v\ annihilation alone can by 
far exceed the observed values in certain cases, they can also be more than one order of 
magnitude below the measurements and the expectations for the secondary spectrum. 

The spectra we just discussed were obtained using the MED set of diffusion parameters 
as given in Table [H together with a solar modulation in the force field approximation with 
4>f = 250 MeV [51]. These diffusion parameters are in fact source of an important uncertainty 
in the flux. For illustration, the plot on the right in Fig. [12] compares the antiproton flux at 
E = O.lm^ obtained with the MED diffusion model (dark blue points) to that obtained with 
the MIN model (light blue points). The scattered points are the allowed scan points of the 
one-sneutrino case, and E = O.lm^ was chosen because as the antiproton spectrum has its 
maximum at an energy that roughly corresponds to 10-15% of the DM mass. The plot gives 
an estimate of the strong dependence of the primary flux on the propagation parameters 
of nearly one order of magnitude just between the MIN and MED model. Similarly, it is 
possible to choose other propagation models (for example the MAX model) that further 
increase the flux by almost an order of magnitude. From Fig. [T2] we see that most of 
our scenarios are below the la uncertainties of the measured flux when assuming the MIN 
propagation model, while strong constraints on our scenarios are expected assuming the 
MED propagation modeljf) 

Let us now turn to the positron fluxes. These are displayed in Fig. [T31 again assuming 
the MIN and MED propagation models and E = m^/10. Here the MIN model gives a 
flux from v\ annihilation that is roughly a factor 1.5 higher than that obtained with the 
MED model. PAMELA has measured the positrons to electron ratio [85] but has not yet 
released the positron fluxes. We can compare however with the secondary flux computed 

3 A detailed fit to the antiprotons and p/p flux ratio will appear elsewhere. 
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Figure 13: Differential ffux for V\ annihilation into positrons at E = O.lm^ including propa- 
gation for two sets of propagation parameters, MIN (light blue points) and MED (dark blue 
points), see Table [TJ 



Model 
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K (kpc 2 /Myr) 


L (kpc) 


V c (km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.7 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 


5 



Table 1: Typical diffusion parameters that are compatible with the B/C analysis [531154]. 

in [86] which lies around <3> e + ~ 5 x 10~ 4 (GeV cm 2 s sr) _1 for the MED propagation model 
in the energy range E = 0.2 — 0.8 GeV relevant for our scenarios; for the MIN model, the 
predictions can increase by roughly a factor 2. From Fig. [13] we see that the positron flux 
from V\ DM annihilation reaches at most $ e + ~ 2 (3) x 10~ 4 (GeV cm 2 ssr) -1 in the MED 
(MIN) model and hence is always smaller than the secondary flux. 

7 Collider signatures of sneutrino DM 
7.1 LHC 

An important characteristic of the light sneutrino DM scenario is the invisible decay of 
the light Higgs [5]. At the LHC, the search for an invisible Higgs will be performed in 
the WW fusion channel with a signature in 2 tagged jets and missing Er- Defining the 
ratio C 2 = &(Hjj)/a(hjj)sM x Br(h —> inv), the region to be probed with £ = 10 fb _1 at 
y/s = 14 TeV corresponds to ( 2 > 0.38 for a Higgs mass below 150 GeV [87]. In our sneutrino 
DM scenarios, the light Higgs is SM-like and the invisible decay h —> v\v\ overwhelmingly 
dominant (« 99%), making for a good Higgs discovery potential in the invisible channel. 

The SUSY signatures also differ from the expectations in the conventional MSSM: while 
squarks and gluinos have the usual cascade decays through charginos and neutralinos, with 
the same branching ratios as in the corresponding MSSM case (see, e.g., [HHIEH]), here the 
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charginos and neutralinos decay further into the v\ LSP. For 90% of the allowed parameter 
points of our scans, both the Xi and the X2 decay to practically 100% into vvi, leading to 
larger missing Et than naively expected for a LSP that weighs only a few GeV. (The Xi 
in fact always decays invisibly.) Note also that in this case q~R — > qx\ and q^ —> qx\ give 
the same signature of jet+E™ lss , differing only in the jet-pr and E™ 1SS distributions. The 
X\ also decays directly into the LSP, with BR(x] t — > / ± i>i) ~ 100% in the large majority 
of the cases. This means that decay chains involving charginos should on average have less 
missing Et than chains involving neutralinos. It also means that X2X1 production leads to 
a single charged lepton rather than the usual trilepton signature. Furthermore if, as we have 
assumed, the LSP is the tau-sneutrino, the charged lepton will be a tau. 

This picture depends only little on most of the parameters that we have fixed in the 
scans, in particular the value of the higgsino mass parameter fi hardly influences the picture 
(lowering \x increases a bit the xf — > W ± Xi decays, but these add to the single lepton 
events, although in a flavour-democratic way). The situation is, however, different for the 
first and second generation slepton masses: lowering them far enough can open up the 
X2 I and/or x\ ~~ >* ^ £r decay channels (£ = e, /x), and this crucially influences the 
experimental signatures. Setting, for instance, 2 = = 100 GeV in the DM allowed 

scan points, while assuming that the RH v e ^ are heavy, gives x° — >* ^ @l decay branching 
ratios of typically up to 30%, and up to 50% if the decay into v T v\ T is suppressed by a small 
mixing angle (the rest goes into LH z/ e;/i 's). Note that decays into staus are mostly absent 
in this case. Taking, as may seem more natural, mj i = m^ 2 « mj 3 , = m^ 2 « 
and |A> e J % |A> T | brings us back to the situation that invisible neutralino decays are 
overwhelmingly dominant. The xf decays democratically into all three lepton flavours in 
this case, xf ~~ >" ^^lh with only a small preference for r ± i>i T . Finally, for A„ — > we 
recover the situation discussed in the previous paragraph, with x® ~^ v tV\ t and Xi ~^ T± ^ir 
having practically 100% branching ratio. 

A detailed study of the LHC potential to resolve the light sneutrino DM scenario, in- 
cluding in particular the determination of the DM mass from q — > q'xf — > q'^^u events, is 
left for future work. 

7.2 ILC signatures 

At an international linear collider (ILC) with y/s = 500 GeV, the main production for the 
light Higgs is e + e~ — > Zh. This allows for precision measurements even if the h decays 
entirely via invisible modes as is the case in our model [90J. The main SUSY production 
processes at the ILC are xfxfi X1X25 X2X21 viP\ as well as f^f^ (and maybe selectrons and 
smuons, depending on the scenario). 

Chargino-pair production has the largest cross-section and should be easily measurable, 
as is illustrated in Fig. HH The peculiarity of the sneutrino DM model is that the chargino 
decay gives a single charged lepton plus missing energy, see the discussion in the previous 
subsection. The cross-section for stau-pair production can reach 70 fb if the staus are kine- 
matically accessible. This occurs, however, only for a small fraction of the successful DM 
models, cf. Fig. [71 The rather heavy staus in our scenarios are a consequence of setting 
m L 3 = m R 3 ■ This assumption was not important for the DM study but has a strong impact 
here. Furthermore, the assumption of heavy selectrons and smuons means that the most 
favourable production processes, notably e + e~ — > e + e~ and e + e~ — > fl + fl~, are kinemati- 
cally not accessible. Relaxing this assumption, for example by assuming universality in the 
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Figure 14: Scatter plot of the unpolarized cross-sections of e + e~ — > XiXi (black) and 7 + inv 
(green) at y/s = 500 GeV as a function of m~±, and of e + e _ — > f^f^ (red) as a function 
of for the allowed scenarios with one light sneutrino. Also shown is the single photon 
cross-section assuming mg R = rri£ L = rrif R . 

slepton masses, could give large selectron/smuon pair-production cross-sections and would 
also affect chargino-pair production, which depends on the mass of the LH v e . 

As mentioned, neutral particles have large branching fractions into invisible states. In 
fact, the x\i x\ ~~ ^ v\v decays both have nearly 100% branching fractions. Therefore, v\v\ and 
XiXj (hj — 1)2) pair production will all contribute to the single photon cross-section. We 
have computed the total single photon cross-section exactly using calcHEP [57] with the cuts 
•p\ > 4 GeV and Obeam-y > 10°. Monte Carlo simulations of this process including backgrounds 
have shown that cross sections around 1.6 fb could be detectable using beam polarisation 
[HUES] B For the allowed points from Section HJ the predictions for the unpolarised cross- 
section are mostly below 1 fb although the cross-section can reach up to 2 pb when M 2 < 
150 GeV, cf. the green points in Fig. HH The dominant contribution arises from the v\v\ 
channel. The mass of the selectron is again a crucial parameter for these processes. To 
illustrate its impact we have also computed the single photon cross-section fixing mg R = 
rrie L = mf L = m fR (instead of 1 TeV). This increases the single photon cross-section by up 
to one order of magnitude, see the blue points in Fig. HH This is due mainly to the increase 
in neutralino production. 

Note finally that although there can be an important additional Higgs production mode 
through the decay i>i — > Dih, this is not detectable since the Higgs decays invisibly. Further- 
more, in our model the v 2 is often too heavy to be pair produced with yfs = 500 GeV. 

8 Conclusions 

In supersymmetric models with Dirac neutrino masses, a weak-scale trilinear term that 
is not proportional to the small neutrino Yukawa couplings can induce a sizable mixing 

4 More detailed analyses including detector simulation for the ILC were also performed for DM masses 
above 100 GeV [551194) . it remains to be seen how can this be applied to our model where particles of different 
masses contribute to the single photon channel. 
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between LH and RH sneutrinos and render the lighter sneutrino mass eigenstate a viable 
dark matter candidate. In particular, the mixed sneutrino can be an excellent candidate for 
light SUSY DM with mass below ~10 GeV, as we have shown in this paper. 

To obtain a very light V\ requires some fine-tuning, as the soft-breaking terms in the 
sneutrino mass matrix are all of the weak scale. Moreover, a viable sneutrino DM candidate 
requires enough mixing to provide sufficient pair-annihilation, while on the other hand the 
mixing should not be too large in order not to exceed the DD limits or contribute too much 
to the Z invisible decay width. In addition, the 1-loop diagrams involving sneutrinos induce 
a negative correction to the light Higgs mass, which can attain a few GeV for large A v . 
In a random scan over 10 5 points for the case of one light sneutrino contributing to the 
DM, we found close to 2000 points that survive all present constraints, including the Z 
invisible decay width, the Higgs and SUSY mass limits, as well as dark matter constraints 
from the relic abundance and the direct detection experiments. These points populate the 
region m Ul ^1 — 8 GeV. We also found that for very light v\ below ca. 5 GeV, small 
M 2 ~ 100 — 150 GeV is preferred in order to enhance V\V\ —> vv annihilation through t- 
channel wino exchange. In the case of three degenerate light sneutrinos, the constrains are 
much stronger and only few points pass all constraints. 

Our results for direct and indirect detection of light sneutrino DM can be summarized 
as follows. First, the cross-section for the spin- independent elastic scattering on nuclei is 
predicted to be <r SI > 10~ 5 pb. That is at most within an order of magnitude of present 
limits for DM masses around 5-10 GeV, and it includes the region favoured by CoGeNT (and 
perhaps CRESST). Second, sneutrino pair annihilation into bb and t + t _ can lead to large 
distortions of the positron and antiproton spectra at low energies. These are being probed 
by PAMELA. Furthermore, the photon flux can be in the range being probed by FermiLAT, 
provided again the zVs annihilate significantly into bb and t + t~. In all cases the signals are 
expected in the low energy range (ca. 0.1 — 3 GeV), a region where measurements are not 
as precise as for higher energies. Third, the neutrino flux resulting from V\ capture in the 
Sun can also be large. However, because the energy of the neutrinos is bounded by the DM 
mass, these scenarios are harder to probe than those with typical weak scale DM. In fact the 
neutrino energy is below the threshold of the large detectors. Only SuperKamiokande has a 
low enough threshold to have some sensitivity to neutrinos from light sneutrino annihilation, 
though these neutrinos do not have enough energy to induce muons that go through the 
detector. One can therefore use only events where muons are stopped in the detector. It is 
also interesting to note that the predicted flux for antineutrinos, and thus for antimuons, is 
often much larger than that for muons. 

While the scenarios with V\ DM that annihilates preferentially into bb have good prospects 
of being detected in the near future, those that have a DM mass below ca. 2 GeV will mostly 
escape detection. Indeed direct detection experiments lack sensitivity for these masses. 
Moreover, the main annihilation channel in this case is into neutrinos, leaving low chances 
for indirect detection in photons, antiprotons and positrons; some possibilities remain if the 
annihilation into t + t~ is significant. Even the neutrino telescopes cannot make use of the 
large branching fraction for annihilation into neutrinos, because the evaporation process in 
the Sun strongly suppresses the neutrino flux for such light UiS. 

Finally a mixed sneutrino LSP leaves distinct signatures in collider experiments. Most 
notably the light Higgs boson and the two lightest neutralinos decay almost exclusively 
into invisible modes, while decays of the lighter chargino give a single charged lepton plus 
missing energy. At the LHC, the typical cascade decays therefore are q — > qxf — > q'l^vi, 
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Ql QX2 ~> Q. v ^i an d Qr — > QXi ~> Qv&i, all giving different amount of missing E T . 
Moreover, X1X2 production gives only a single charged lepton. At the ILC, z>i£* and 
= 1,2) pair production will all contribute to the single photon cross- sect ion. 

The details of the DM signatures at colliders depend of course on the assumptions that 
are made on the rest of the spectrum, in particular the first and second generation of slep- 
tons. The predictions in astroparticle experiments have uncertainties from astrophysical and 
nuclear parameters. These include uncertainties in the quark content of the nucleon, on the 
local dark matter density and velocity distribution, on the dark matter halo profile and, for 
charged cosmic rays, on the parameters of the propagation model. 

Last but not least if signals are found in astroparticle and collider experiments, the 
challenge will be to determine the precise DM properties and the underlying new physics. 
This requires in particular collider measurements of the masses and couplings of the DM and 
other new particles associated with it in order to refine and test the theoretical predictions 
for astroparticle observables. Moreover, if an E™ 1SS signal is seen at the LHC, it should be 
confirmed in direct DM detection and vice versa, the agreement of the DM mass and cross- 
section determined in the two ways being a crucial test. To this end it will be interesting 
to investigate how well our light sneutrino DM scenarios can be resolved at the LHC by 
exploiting the q — > qxf — > q'^vi cascade decay. Such a signature has been suggested in [95] 
as a method of measuring mass differences in the MSSM. This will be a future work. 
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A Vertices involving sneutrinos 

The Feynman rules of relevant V\ couplings include 
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Here is the charged lepton mixing matrix that is diagonal in flavour space, N is the 
neutralino mixing matrix and U, V the chargino mixing matrices. We use the SLHA notation 
for the MSSM part of the Lagrangian 
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